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Precise lattice parameter measurements in single crystals are achievable, in principle, by X-ray 
multiple diffraction (MD) experiments. Tiny sample misalignments can compromise systematic 
usage of MD in studies where accuracy is an important issue. In this work, theoretical treatment 
and experimental methods for correcting residual misalignment errors are presented and applied to 
probe the induced strain of buried InAs quantum dots on GaAs (001) substrates. 
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I. INTRODUCTION 

Opto-electronic devices based on self-organized InAs 
quantum dots (QDs) are very suitable for metropolitan- 
area network applications. Their wavelength emission 
range match the absorption minima of the optical fibers; 
this optical property is still under improvement to pro- 
vide an even better emission- absorption match. 1 The ac- 
tual challenge resides in increasing the optical efficiency 
of these devices, which is possible by increasing the den- 
sity of optically active QDs. Recent studies have demon- 
strated drastic reduction in the number of optically active 
structures after the growth of a few atomic GaAs lay- 
ers over the QDs. Since surface probe techniques are no 
longer useful to inspect the physical structure of buried 
QDs; alternative structural characterization procedures 
have become relevant. 

X-ray Renninger scanning 5 (XRS) is the most precise 
technique for accurate lattice parameter determination 
in single crystals.- However, in practice, even in very 
good equipments and careful alignment procedures, some 
tiny sample misalignments are always present. They 
have compromised XRS application in studies where ac- 
curacy is an important issue. In this work, we describe 
how to consistently correct systematic errors in XRS to 
achieved high sensitivity in the lattice-parameter varia- 
tion (Aa/a 10 -5 ). Combining such an accuracy with 
the shallow penetration depth of the X-ray wavefield un- 
der Bragg-surface diffraction^ — a particular multiple 
wave diffraction case — the residual average strain field 
on the substrate lattice of devices with buried QDs has 
been characterized along two orthogonal in-plane direc- 
tions. 



II. X-RAY DIFFRACTION GEOMETRY 

X-ray multiple diffraction (MD) in crystals are excited 
when the incident beam direction, wavevector fc, fulfill 
two conditions summarized by 




FIG. 1: Three-beam X-ray diffraction in crystals with the 
primary diffraction vector P aligned along the scanning axis 
0, P || z. When a secondary reflection S diffracts, C = P — 
S provides the energy coupling S P reflection. Rotation 
matrices R X (S X ) and R y (5 y ), account for tiny misalignments 
of diffraction vectors with respect to the (x, y, z) goniometer's 
reference system, k is the incident beam wavevector. 



and 



fc-S = -SS/2. 



(2) 



k P = -P P/2 



(1) 



P and S are the diffraction vectors of the primary and 
secondary reflections, respectively. In XRS, the primary 
reflection is kept excited, i.e. Eq. (QJ fulfilled, while the 
crystal rotates around P. It requires P parallel to the 
goniometer cf) axis, 4 the z axis in Fig. 1. MD occurs when 
secondary reflections are excited by the <p rotation and, 
therefore, Eq. @ is also fulfilled. The primary intensity 
changes as secondary reflections diffract, e.g. Fig. 2. 

The most well known expressions for determining MD 
positions in XRS were obtained from Eq. (|2j! 6,T assuming 
that the primary reflection is always aligned, i.e. Eq. (QJ 
fulfilled during a complete <j) rotation of 360°. An alter- 
native approach to determine such positions without the 
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P || z constrain has been obtained 8 by writing P = S + C 
in Eq. which leads to 

k • C = -C • C/2 - C • S. (3) 

Sample misalignments S x and S y , as shown in Fig. 1, 
generate small Aco = uo — uoq and A</> = <\> — fa corrections 
in the incident beam direction k(uj : fa) ~ fco + Auk u + 
Afikcf). These corrections can be determined by a trivial 
system of linear equations 



ku • S k^ • S 
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(S/2 + fco) • S 


k^j • C k^ • C 








(C/2 + S + fco) -C 



(4) 



derived from Eqs. (0) and (0). k u = dk/duo and k<j> = 
dk/d(j) are calculated at c<;o and </>o, the incidence and 
azimuthal angles for exciting the MD in non-misaligned 
samples, i.e. at fco = fc(^o, fa)- Moreover, G = Go + AG 
stands for the S or C diffraction vectors where 

AG = (G z 5 y , G Z 5 X , — G x 5 y — G y S x ) 

and G = (G x , G y , G z ). 

In non-misaligned sample the MD positions are ob- 
tained from Eq. alone, which provides 

\ a A|S|/2 - S^sincjo 
cos (0o — a) = cos jj = (5) 

S^COSCJo 

where S x = S xy cosa = x • So, S y = S^sina = y • So, 
and uq is the Bragg angle of reflection P. 

There are two azimuthal positions where the same sec- 
ondary reflection is excited: fa = a — (3 and fa = a + /3. 
Lattice-parameter determination is based on the exper- 
imental measurements of both positions. They provide 
2/3 e xp — fa,exp — 4>i,exp for the calculation of the unit cell 
parameters from Eq. (0, as described elsewhere (see for 
instance Ref . 0) . The residual sample misalignments are 
measured by rocking curves of the primary reflection at 
= 0, 90°, 180°, and 270°. If cj , ^ 90 , ^iso and co 2 70 
are the respective rocking-curve peak positions after the 
final alignment, 

S x = (^270 - ^9o)/2 and S y = (c^i 8 o - ^ )/2. (6) 

The A0 n= i 5 2 corrections in the azimuthal positions are 
obtained from Eq. Q, and they can be used either to 
refine the 2j3 exp value according to 

^ftexp = (fa,exp A0 2 ) (fa,exp ^fa) (7) 

or to estimate misalignment effects on different secondary 
reflections, for instance those providing some useful in- 
formation about the crystalline structure of the sample. 
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FIG. 2: XRS of the 002 GaAs reflection carried out with 
X-ray wavelength A = 1.330234A. The peaks at <f> = 175.08° 
(No. 4) and 184.92° (No. 5) occur when the 111 and 111 
secondary reflections are excited, respectively. MD positions 
(closed circles, 1 to 8) owing to other secondary reflections of 
the 111 family are indicated in the inset; their measured <j> 
positions are provided in Table 1. 



III. EXPERIMENTAL 

Data collection has been carried out at the Brazilian 
Synchrotron Light Laboratory (LNLS) D12A (XRD-1) 
beam-line, with the polarimeter-like diffractomer; 4 small- 
est step size is 0.0004° in both uj and <\> axes. The wave- 
length, A = 1.330234A from a Si (111) double-crystal 
monochromator, was measured by rocking-curves of the 
111 and 333 silicon reflections. Vertical incidence plane 
scattering upwards (a-polarization). The sample is a 
commercial GaAs (001) substrate with In As quantum 
dots (QD) grown by molecular beam epitaxy (MBE) us- 
ing a rate of 0.009 monolayer per second (mL/s). There 
is a 300A thick GaAs cap layer. 



IV. RESULTS AND DISCUSSIONS 

Tabled shows the fa and fa positions of the 111, 111, 
111, and ill secondary reflections in the XRS of the 002 
GaAs reflection. The peak positions were determined 
by fitting the intensity data with lorentzian-gaussian 
convolution curves: a lorentzian function standing for 
the intrinsic profile (FWHM ~ 0.0048°) while a gaus- 
sian accounts for the instrumental broadening (FWHM 
~ 0.0060°). A short portion of the XRS around = 180° 
is shown in Fig. 2; the [110] direction is the reference for 
= 0. To check for instrumental errors, 0-scans of the 
peaks 4 and 5 were repeated several times after ±360° 
rotations in fa the observed differences were not larger 
than 0.00075°. Each <j> scan has been performed at the 
maximum of the (002) rocking-curve, which was carried 
out about 0.5° before the MD peak positions. 
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FIG. 3: Experimental (open circles) and simulated (solid 
lines) rocking curves, o;-scans, of the 002 primary reflection 
measured at = 0, 90°, 180° and 270° (top- left of each 
scan). The scans were simulated by two gaussian functions 
(dashed lines) whose splitting is about Acu ~ -0.00422(24)°, 
and hence Ac/c ~ 3.0(2) x 10" 4 . 6' x (= 0.0024°) and 
Sy (= 0.0103°) are the expected misorientations according 
to the P values in Table [I] 



The rocking curves used to characterize the resid- 
ual misalignments are shown in Fig. 3; they provide 
S x ~ 0.0009° and S y ~ 0.0067°. Since these values 
are very small, Eq. Q can be linearized by numeri- 
cal derivation, and the azimuthal corrections written as 
A(j) n = A n 5 x + B n S y . The 4-fold axis symmetry of the 
measured secondary reflections establish some relation- 
ships among the A n and B n coefficients, as given in Ta- 
ble [I] (last two rows). By replacing them into Eq. {7J, it 
is possible to demonstrate that 
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FIG. 4: Schematic illustration of the compressive stress gen- 
erated by InAs QDs on GaAs substrate. There is also a 
300A thick GaAs cap layer over the QDs (not shown in the 
figure) . 



Note that |S| = |So| is invariant under any rotation, but 
their components are not. Therefore, S x , S y and S z in 
Eq. © must be calculated for x = (1, 1, 0)/<\/2, y = 
(1, 1, 0)/a/2 andz = (0, 0, 1), i.e., S x = (ft + k)/(ay/2), 
S y = (-ft + k)/(ay/2) and S z = 1/c. 

The observed tetragonal deformation v ~ 1.77 x 10 -4 , 
corresponds to an average value in a shallow layer just 
below the surface, not thicker than 0.3 fim^ where the 
induced strain due to the InAs QD is significant. Al- 
though the InAs epitaxial growth generates an expansive 
stress in the substrate lattice under the QDs, the adja- 
cent regions are compressed as schematically illustrated 
in Fig.QJ This in-plane compressive strain may have been 
propagated into the epitaxial cap layer. 

An error in the wavelength value of AA = +0.0015A 
could be responsible for the j3 values obtained in Table 1; 
but, the used method for determining A assures a preci- 
sion that is at least 10 times better, AA/A < 2.0 x 10 -4 . 

Differences in the v values along the [110] and [110] 
directions, i.e. Alio] ^ Alio]' would yield a unit cell 
slightly twisted near the substrate surface. The variation 
in the 7 angle, between the a and b lattice vectors, is 
estimated from Eq. (0 as 



Ano] =09iii+#EIi)/2 (8a) 

and 

Alio] = (/?in+ An )/2 (8b) 

do not depend on S x and S y . In other words, the aver- 
age /3 exp for the secondary reflections with in-plane com- 
ponents along the [110] and [110] orthogonal directions, 
/%io] an d Alio] , are misalignment-free experimental val- 
ues and, therefore, useful for precise lattice parameter 
measurements in both directions. 

Assuming a tetragonal substrate lattice distortion, its 
in-plane lattice parameter is given by a = b = (1 — u)ao 
while, to preserve the unit cell volume, c = (1 + 2v)clq. 
For a = 5.6534A, |P| = 2/c, S = (ft/a, fe/6, 1/c), and 
ft, k = ±1, Eq. J5J) provides 

(3{v) = 85.08984° - 33.85°^ - 5.2°^. (9) 



Alio]" Ano] = (Clio -Cno)A 7 (10) 

where C I10 = 0.1276 and C no = -0.1302; then 7 = 
89.9985°±0.0018°. Non-uniform self-organization of QDs 
plus the movement of the beam spot on the sample sur- 
face during the XRS could produce similar results. It 
could also be responsible by the observed disagreement 
between the S XyV and 5' x values. 



V. CONCLUSIONS 

In summary, refinement of XRS data for residual sam- 
ple misalignment errors has led to a precise tool for prob- 
ing in-plane lattice strains, as small as Aa/a — 10 -5 . Al- 
though, the strain of InAs QDs on the surface of GaAs 
(001) substrates was the subject investigated here, the 
method can be extended for studying other epitaxial 
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TABLE I: Azimuthal <j> positions of the 111, 111, 111, and 111 secondary reflections in the XRS of the 002 GaAs reflection. 



Each position was measured three times (rows 1, 2 and 3) as explained in the text, 



: + <fimin)/2, £ = (4>n 



(fimin ) /2 5 



(02 — 0i)/2, /3 are the average misalignment-free values, v stands for the unit-cell tetragonal distortion [Eq. @] and 



Pe 

p m%s =/3+ (A0 2 - A0i)/2 [Eq. Q] where A0 n = A n 6' x + B n 6' y , S' x = 0.0024°, 
by numerical derivation of Eq. Q. Angular values are given in degrees. 



and 8' — 0.0103°. A n and B n were estimated 



111 



111 



111 



111 



1 

2 
3 

4> 
e 

ftexp 

P 

V 

A n 

B n 



-85.07750 
-85.07716 
-85.07683 
-85.077165 
±0.000335 



85.08599 
85.08575 
85.08548 
85.085735 
±0.000255 



94.91867 
94.91853 
94.91819 
94.91843 
±0.00024 



265.09121 
265.09141 
265.09209 
265.09165 
±0.00044 



85.081450±0.000295 



85.08661±0.00034 



85.08403±0.00032 
(1.716 ±0.096) x 10~ 4 
85.0815 85.0865 
0.0207 0.0207 -0.0207 -0.0207 

0.2412 -0.2412 -0.2413 0.2413 



4.92085 
4.92102 
4.92105 
4.92095 
±0.00010 



175.08692 
175.08688 
175.08730 
175.08709 
±0.00021 



184.92072 
184.92061 
184.92064 
184.920665 
±0.000055 



355.08935 
355.08881 
355.08880 
355.089075 
±0.000275 



85.083070±0.000155 



85.084205±0.000165 



85.08364±0.00016 
(1.832 ±0.048) x 10~ 4 
85.0831 85.0842 
0.2412 -0.2412 -0.2413 0.2413 

-0.0207 -0.0207 0.0207 0.0207 



nanostructures by their strain field on the substrate lat- 
tice. 
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